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ABSTRACT 

We investigate the evolution in the dark matter halo circular velocity function, measured at a fixed physical 
radius of 20 kpc (V20), which is likely to be a good proxy for galaxy circular velocity, in the Millennium-II 
simulation. We find that the V20 function evolves remarkably little since z ps 4. We analyze the histories of the 
main progenitors of halos, and we find that the dark matter distribution within the central 20 kpc of massive 
halos has been in place since early times. This provides evidence for the inside-out growth of haloes. The 
constancy of the central circular velocity of halos may offer a natural explanation for the observational finding 
that the galaxy circular velocity is an excellent predictor of various galaxy properties. Our results also indicate 
that we can expect a significant number of galaxies with high circular velocities already at z = 4 (more than 
one per 10 6 /i _3 Mpc 3 with circular velocities in excess of 450 km s _1 , and more than one per 10 45 /i _3 Mpc 3 
with circular velocities in excess of 350 km s -1 ). Finally, adding baryonic mass and using a simple model 
for halo adiabatic contraction, we find remarkable agreement with the velocity dispersion functions inferred 
observationally by Bezanson et al. (201 1) up to z pa 1 and down to about 220 km s _1 . 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: formation — galax- 
ies: abundances — galaxies: halos 



1. INTRODUCTION 

Simulations of individual dark matter halos have provided 
evidence for the inside-out growth of halos in a ACDM cos- 
mological model. These studies have shown that the inner 
few tens of kpc remain almost unchanged since high redshift 
both in massive halos (e.g. Loeb & Peebles 2003; Gao et al. 
2004) and in halos expected to host Milky-Way sized galax- 
ies (Wechsler et al. 2002; Diemand et al. 2007; Wang et 
al. 2011). From studying the mass-concentration relation of 
halos in cosmological simulations, it has become clear that 
the central density of haloes is correlated with their epoch of 
formation, which provides additional evidence for inside-out 
growth (Bullock et al. 2001a; Wechsler et al. 2002; Zhao et 
al. 2003). 

The almost constant central densities of halos may have di- 
rect consequences for galaxy formation: At z = 0, the circu- 
lar velocity at 10 kpc for dark matter halos has been found to 
be a good proxy for galaxy circular velocity and galaxy ve- 
locity dispersion (Loeb & Peebles 2003; Trujillo-Gomez et 
al. 2011). It is likely that this is also the case at higher red- 
shift, and that central densities of halos can be used to predict 
galaxy circular velocities. 

Motivated by these findings, we study the evolution of the 
circular velocity of dark matter halos at 20 kpc (hereafter i^o) 
in this work. We choose a radius of 20 kpc because galaxies 
in the nearby Universe are thought to have nearly isothermal 
mass profiles out to this scale. The scale is also small enough 
that we can expect structure to have formed very early (e.g. 
Diemand et al. 2007), and at the same time large enough that 
it is not completely dominated by baryonic processes. 

We extend previous studies, which have focussed on in- 
dividual halos, by studying the same quantity for a statisti- 



cal sample of dark matter halos in the Millennium-II high- 
resolution cosmological simulation. Thus, we can for the first 
time explore the evolution of the small-scale circular velocity 
function of halos as a function of redshift. 

2. DATA 

We use the publicly available dark matter data from the 
Millennium-II simulation (hereafter MS -II, Boylan-Kolchin 
et al. 2009), which is contained in a box of side 137 Mpc 
and has a particle mass of 9.45 • 1O 6 M . According to eq. 20 
in Power et al. (2003) the MS-II mass resolution is sufficient 
to resolve the circular velocity profile of halos down to 20 kpc 
and i? max for halos with masses above 10 10 M©. Particle data 
for this simulation is not available. Thus, we calculate V20 
from -R m ax, -R200 an d W200 analytically, assuming a Navarro, 
Frenk & White (1997, NFW) profile: 



p(R) = 



Po 



(R/R s )(l + R/R s 



(1) 



with the scale radius R s and po as the free parameters of the 
profile. The concentration of a NFW halo is defined as the 
ratio between scale radius and virial radius, c = i? s /i?2oo- 
Given that the maximum circular velocity is reached at a ra- 
dius i? max = 2.163 • R s for a NFW halo, it follows that 

c= 2.163 -i?2ooAR max ■ (2) 

Integrating the NFW profile, the enclosed mass at radius R is 
given by 



M (R) = AnpoRs 



ln(l + cx) 



cx 
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cx 



(3) 
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with x = R/R 200 - Given that v cilc {R) = \J GM( -< R ^ 
then calculate the circular velocity at 20 kpc using 

1 ln(l + cx) 



V20 = V200 
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FIG. 1 . — Evolution of the density profile of a randomly selected halo with a (2 = 0) mass of 7 ■ 10 12 Mq, and with a virial radius of 370 kpc (left) and with a 
mass of 2.8 ■ 10 13 Mq and 610 kpc (middle panel). Crosses show the location of the virial radius at each redshift. In the right panel, we show the average profile 
of 22 halos in a narrow bin in halo mass from 6.8 ■ 1O 12 M0 to 7.5 ■ 10 12 Mq, following the main progenitor history of halos. This figure shows that the density 
profile within « 20-50 kpc undergoes only little evolution since 2^4. 

with results from the Milli-MSI simulation. In Fig. [2] we 
show the halo mass, virial velocity, maximum circular veloc- 
ity («max) and V20 distribution for both simulations, with V20 
calculated directly from the particle data for the Milli-MSI, 
and indirectly for the MS -II, at z — and z = 2. In principle, 
the inner profile of halos can only be expected to be well re- 
solved down to the limit given by Power et al. (2003). If we 
demand that less than 10 % of Milli-MSI halos fall below this 
limit, the resolution limit is about w V j r =250 km s" 1 at z = 0, 
and w V i r =117 km s _1 at z = 2. We find that V20 for the two 
simulations is in reasonable agreement, although it appears to 
be slightly lower in the Milli-MSI simulation. 

3.2. Evolution in V20 

In Fig. [3]we show the evolution of the halo mass function, 
virial velocity function, v max function and V20 function. The 
halo mass function evolves dramatically over this time period. 
The virial velocity function and the v max function also show 
significant evolution at the high mass end, indicating the for- 
mation of groups and clusters of galaxies. At intermediate 
values of w max « 200 km/s, little evolution is seen, in good 
agreement with results by Bullock et al. (2001b). Finally, 
very little evolution in the entire range of velocities is seen for 
the V20 function, which is well in place since z = 4. As we 
show below, this is likely connected to the inside-out growth 
of halos. 

We then consider the main progenitor histories of halos in 
narrow halo mass bins at z = in Fig. |4j averaging over 
100^-00 central halos in the three lower bins, and over 43 
halos in the highest mass bin. In agreement with the previous 
figure, we can see circular velocities growing more slowly 
than halo mass, particularly at late times. Since 2 = 1, neither 
the average V20 nor the average t> V i r has changed significantly 
in all mass bins. For the three higher bins, V20 has even stayed 
practically constant since z = 2. This figure is interesting 
since it allows us to directly read off the expected value of V20 
for the main progenitor of a halo with a given final mass at 
z = 0. A halo with a mass of 1O 13 M , which is expected to 
host a galaxy with a stellar mass of about lO n M today (e.g. 
Moster et al. 2012) will have had V20 ~ 250 km s _1 since 
z = 2, and had reached V20 ~ 160 km s _1 already at z = 6. 
A halo with a mass of 1O 12 M , which is expected to host a 



with x=2Q kpcAR vil .. Klypin et al. (2001) show that NFW 
profiles fit numerically simulated profiles with a relative error 
of below about 10 % for radii larger than 1 % of the virial 
radius. The density at 20 kpc should thus be well described 
by a NFW profile. 

For most of our analysis (with the exception of Fig. [T]and 
2] where we only consider centrals), we use both central and 
type 1 satellite galaxies, only leaving out the 'orphan' galaxies 
without associated subhalos. For the satellite galaxies, we use 
-Rmax and 1*200 at the time of infall to calculate V20, i-e. we 
assume that V20 stays unchanged after infall. This assumption 
seems reasonable, see Diemand et al. (2007). 

We also make use of the lower-resolution Milli-Millennium 
simulation (hereafter Milli-MSI, Springel et al. 2005), which 
has a particle mass of 1.18 • 10 9 M Q . For this simulation, in- 
dividual particle data are available, which we use to calculate 
individual halo profiles, and to compare with the results from 
the MS-II simulation. Both the MS and MS-II simulations are 
based on a WMAP1 cosmology (Spergel et al. 2003). 

3. RESULTS 
3.1. Particle data 

As a first step, we consider the dark matter profiles of cen- 
tral halos in the Milli-MSI simulation, calculated from parti- 
cle data. In Fig. [T] left and middle panel, we show the evo- 
lution of the density profile up to z — 4 for two randomly 
selected halos. The halo shown in the left hand panel reaches 
a virial mass of 2.8 • 1O 13 M and a virial radius of 610 kpc, 
the one in the middle panel a mass of 7 • 1O 12 M and a virial 
radius of 370 kpc by z — 0. In the right panel of the same fig- 
ure, we show the average profile of 22 halos in a narrow bin in 
stellar mass from 6.8 • 1O 12 M to 7.5 • 1O 12 M . The central 
parts of haloes (inside 20-50 kpc) show only relatively little 
change, while there is significant growth on larger scales. We 
also show the location of the virial radius and virial density 
at different redshifts as crosses to indicate the significant evo- 
lution of those quantities. As a side note, absence of change 
in the inner part of the profile does not mean that the same 
dark matter particles stay in the inner region of the halo since 
z « 4 (e.g. Gao et al. 2004). 

As a next step, we check the agreement between the sta- 
tistical distribution of halo properties in the MS-II simulation 
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FIG. 2. — We compare results from MS-II (lines) to results from the Milli-MSI simulation particle data (symbols) down to the resolution limit. Filled symbols 
are results above the resolution limit according to Power et al. (2003), empty symbols below. The similarity of the results shows that velocity profiles of halos 
are reasonably well converged even below the resolution limit by Power et al. (2003). 

stellar mass as given in Moster et al. (2012). For the cold gas 
fraction, / gas = M gas /M», we use the average f gas (z,M*) 
from the Guo et al. (201 1) semi-analytical model. The result 
is shown in the top middle panel of Fig. [5] The difference to 
the dark matter-only model is small. 

As a next step, we take into account the effect of adiabatic 
contraction (AC) using the Blumenthal et al. (1986) formula, 
iteratively solving the equation 



galaxy with a stellar mass of about 2.5 • 10 M Q at z = 
will have had V20 ~ 180 km s _1 since z = 2, and W20 ~ 
100 km s _1 already at z = 6. From this we conclude that 
the main progenitors of today's massive galaxies will likely 
exhibit high circular velocity already at high redshifts. 

3.3. Comparison to observational results 

As a next step, we compare our findings to observational 
results of the galaxy velocity dispersion a by Bezanson et 
al. (2011) and Bernardi et al. (2010), assuming 0= t | 2o/v / 2- 
Bernardi et al. (2010) calculate velocity dispersion functions 
from the SDSS using the spectroscopically measured quantity 
within the SDSS fiber. Bezanson et al. (201 1) estimate the ve- 
locity dispersion of higher redshift galaxies using the stellar 
masses and radii, using a method calibrated onz = data. 

Assuming no halo contraction and no baryonic content, the 
agreement between the MS-II results and the observations is 
already encouraging (Fig. [5] top left panel). But of course, 
the presence of baryonic matter in the central regions of the 
halo, and the potential effect of halo contraction (Blumenthal 
et al. 1986; Trujillo-Gomez et al. 201 1) should both serve to 
increase the central circular velocity dispersion. 

We add baryonic mass according to 

M 20 kpc = M20kpc,DM • (1 - /b) + M* • (1 + /gas) (5) 

with /b = 0.17 the cosmic baryon fraction. We estimate M* 
from the redshift-dependent relation between halo mass and 



n ■ M tot ( n ) = n ■ [M tot ( n ) ■ (1 - f b ) + M* • (1 + / gas )] (6) 

for n, with rf =20 kpc. In case this formula results in halo 
expansion instead of contraction, we assume that the halo pro- 
file simply stays fixed. While the magnitude of the effect of 
adiabatic contraction, and even its existence, is still unclear 
and might differ for different galaxy types (e.g. Gnedin et al. 
2004; Dutton et al. 2011), taking into account both the con- 
tribution of baryonic material and AC results in near perfect 
agreement with the Bezanson result (top right panel of Fig. 
|5]l at the high velocity dispersion end. At lower velocity dis- 
persions, the agreement breaks down. This is similar to the 
results obtained by Trujillo-Gomez et al. (201 1) at z = 0. 

The disagreement between our predictions and the results 
of Bezanson et al. (2011) at lower velocity dispersions may 
indicate that the simple relation between velocity dispersion 
and circular velocity breaks down for disky galaxies with a 
small bulge component. In addition, the estimates for veloc- 
ity dispersions used by Bezanson et al. (201 1) might become 
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FIG. 3. — Evolution of the halo mass, virial velocity, n ma x and V20 functions at different redshifts. While the halo mass function undergoes substantial growth 
at all masses, the t> ma x and virial velocity function only grow at the high velocity end, and actually decrease at the low velocity end. The V20 function undergoes 
remarkably little change since 2 = 4. 



progressively uncertain for lower-mass galaxies, due to the in- 
creasing contribution of face-on disks and the higher gas and 
dark matter content of low-mass galaxies (but see Taylor et 
al. 2010 who find that the method works well down to 100 
km s _1 at z = 0). The fact that the agreement improves to- 
wards lower redshifts might indicate that velocity dispersions 
are particularly difficult to estimate for highly star-forming, 
gas-rich galaxies. Bezanson et al. (2012) show that the a be- 
low which more than 50 % of all galaxies are star-forming 
decreases towards lower redshift, a signature of downsizing. 
In all, the cause of the disagreement at low velocity disper- 
sions is not entirely clear, but the agreement is excellent in 
the high-mass regime where elliptical, passive and probably 
gas-free galaxies dominate. We note that previous studies 
have found strong and puzzling disagreement between obser- 
vations and model predictions at even lower velocity disper- 
sions, below 100 kins" 1 (Zavala etal. 2009), with many more 
low-velocity dispersion galaxies predicted by models. One 
potential explanation is that dark matter halos are strongly in- 
fluenced by baryonic processes in this regime (Sawala et al. 
2012). 

Finally, in the bottom row of Fig. [5] we show the cumula- 
tive distributions of V20 and w max . In the bottom left panel, 
we show calculated from dark matter alone, in the mid- 
dle panel v^o including the effect of baryons and AC. These 
two figures can be used to read off predictions about the num- 
ber density of galaxies with V20 above a certain limit. In the 



bottom right panel of the same figure, we compare our cumu- 
lative t> max function to the analytic fit of Klypin et al. (201 1) 
based on the Bolshoi simulation, which uses a WMAP7 cos- 
mology (Komatsu et al. 2011). The agreement is reasonable 
given the different cosmology, with the MS-II simulation pre- 
dicting more high mass haloes, as expected. Similar results 
have been obtained from a comparison between WMAP1 and 
WMAP7 velocity functions in Guo et al. (2013). 

4. SUMMARY AND DISCUSSION 

Using the MS-II simulation, and particle data from the 
Milli-MSI, we have investigated the evolution of the circular 
velocity of dark matter halos at 20 kpc, V2o- Our main results 
are the following: 

• For halos with mass log(M/M Q ) > 12 , the mass pro- 
file within 20-50 kpc changes only little since about 
z « 4, despite significant growth in halo mass. Fol- 
lowing the main progenitor history, we find almost no 
evolution in V20 over this time period. 

• As a consequence, the U20 function of halos is almost 
in place at z = 4. 

• By using V20 as a proxy for the circular velocity of 
the galaxy, we can reproduce observations of the ve- 
locity dispersion of massive galaxies by Bezanson et al. 
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FIG. 4. — Evolution of the average halo mass, virial velocity and Umaxi an d v 20, for the main progenitors of halos in four bins of M200 at 2 = 0. The top two 
panels and the bottom right are based on dark matter results only, the bottom left panel shows V20 after correction for baryons and AC. Results for V20 are only 
shown as long as V20 < f max ■ The right axis of the bottom two panels show rescaled results for a = V20/ V2. Confirming results from Fig. [3] V20 undergoes 
only little change since 2 ss 2 — 4, depending on the halo mass bin. 



(2011) remarkably well. The agreement is nearly per- 
fect if we take into account the baryon content of halos 
inferred from abundance matching, and by considering 
the Blumenthal et al. (1986) prescription for adiabatic 
contraction. The observations of Bezanson et al. (2011) 
are thus in excellent agreement with basic ACDM pre- 
dictions on the evolution of structure in the Universe. 

These results have interesting implications. First of all, we 
expect a significant number of galaxies (more than one per 
10 6 ft~ 3 Mpc 3 ) with circular velocities in excess of 450 km 
s _1 already at z = 4. We also expect more than one per 
10 45 /i _3 Mpc 3 with circular velocities in excess of 350 km 
s _1 up to z = 4. Our results thus predict that high line widths 
should be found in high redshift galaxies with the Atacama 
Large Milimeter/Submilimeter Array (ALMA). 

Second, we infer that the main progenitors of today's most 
massive galaxies must have had high circular velocities al- 
ready at high redshifts. This means that the circular velocity 
can be used to link massive galaxies with their progenitors and 
descendants. It also helps to understand the fact that high red- 
shift galaxies have higher circular velocity at given dynamical 
mass than their local counterparts. Recently, van de Sande et 
al. (201 1) have measured a velocity dispersion of ps 300 km 
s _1 in a galaxy at z = 1.8 with M^/Mq w 1.7 • 10 11 , corre- 
sponding to a very high velocity dispersion for its dynamical 
mass compared to today's galaxies. As seen in Fig. [4] the 
galaxy will continue to grow in halo mass (and probably also 
in stellar mass) to z = 0, but its circular velocity has likely al- 
most reached its final value by z w 2. An even more extreme 
example is the extremely high velocity-disperson galaxy (510 



km s _1 ) discovered by van Dokkum et al. (2009). 

Finally, our result might offer an explanation for a puzzling 
recent observational result. Wake et al. (2012a, b) have found 
that galaxy circular velocities are an excellent predictor not 
only of the clustering properties of galaxies, but also of galaxy 
colours, specific star formation rates and ages, being superior 
to galaxy stellar mass or surface mass density. As we have 
shown in this work, the galaxy circular velocity, unlike the 
galaxy stellar mass, stays almost constant since early times. 
It thus might constitute a galaxy property that is more funda- 
mental than stellar mass, retaining some memory of the initial 
conditions of the galaxy's formation. 

We thank the referee for helpful comments, which 
improved the paper. We acknowledge funding from 
ERC grant HIGHZ no. 227749. SQL databases 
containing the simulation data from Millennium and 
Millennium-II simulations are publicly released at 
http : / / www . mpa-garching . mpg . de /millennium. 
The Millennium site was created as part of the activities of 
the German Astrophysical Virtual Observatory. Tables for 
the velocity function of halos as shown in this paper can 
be found at http://www.strw.leidenuniv.nl/ 
galaxyevolution/ v2 function/. 
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FIG. 5. — Top panels: Comparison of a at 20 kpc (i.e. t'2o/v / 2) to results by Bezanson et al. at z = 0.45 and z = 1.05, and Bernardi et al. (2010) at z = 
(crosses). We also show predictions for z = 2 and z = 3 as dashed lines. Results are shown for dark matter alone (top left panel), dark matter and baryons (top 
middle panel), and dark matter, baryons and AC (top right panel). Agreement between theoretical predictions and the Bezanson et al. (201 1) results are excellent 
for the high mass end of the velocity dispersion function. Bottom panels: Cumulative functions. Bottom left panel: ^20-function, only taking into account dark 
matter. Middle panel: Same like left hand panel, but with the contribution of baryons and AC added. Right hand panel: Cumulative « max -function, compared to 
the analytical fits of Klypin et al. (201 1) based on the Bolshoi simulation. 
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